Crystal-field effects in the mixed-valence compounds Yb 2 M 3 Ga9 (M= Rh, Ir) 
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Magnetic susceptibility, heat capacity, and electrical resistivity measurements have been carried 
out on single crystals of the intermediate valence compounds Yb2Rh3Gag and Yb2lr3Gag. These 
measurements reveal a large anisotropy due apparently to an interplay between crystalline electric 
field (CEF) and Kondo effects. The temperature dependence of magnetic susceptibility can be 
modelled using the Anderson impurity model including CEF within an approach based on the Non- 
Crossing Approximation. 

PACS numbers: 75.30.Mb,75.20.Hr,71.27.+a,71.28.+d 



The intermediate valence compounds pose one of the 
most challenging problems of strongly correlated elec- 
tron systems. Different ingredients contribute to the 
complexity of these fascinating systems: the presence 
of strong Kondo interactions, the level structure of the 
crystal electric field (CEF) /-orbitals, the different hy- 
bridizations between each level and the conduction band, 
and the eventual coherence effects and magnetic interac- 
tions introduced by the periodicity of the Kondo latticed 
Strong valence fluctuations are observed in the inter- 
metallic compounds with Ce, Yb and U. In particular, 
Yb compounds attract a great deal of interest because 
the trivalent Yb ion is at least in some sense the hole 
counterpart of the Ce 3+ ion which has one electron in its 
4/ shell. As in the case of the Ce compounds^ the Yb- 
based intermetallics exhibit a diversity of physical prop- 
erties that remain to be understood^ 

The isostructural series of compounds R2M3X9 (R = 
La, Ce, Yb, U; M=Co, Rh, Ir; X = Al, Ga) exhibit an- 
tiferromagnetic ordering for R=Yb, X=A1, and mixed- 
valence behavior for R=Ce,Yb and X=Gai£&! All the 
U-based compounds order antiferromagnetically at tem- 
peratures below 40 Ki4 The Yb2MsGag compounds are 
a suitable class of materials for studying the difference 
in the electronic structure between the magnetically or- 
dered Kondo lattice and the mixed-valence systems. The 
orientation-dependent temperature T max of the maxi- 
mum in the magnetic susceptibility suggests the possibil- 
ity of an anisotropic Kondo effect &2i Previously, Petro- 
vic et al. studied ternary R-Ir-Ga compounds (R=rare 
earth) that were assigned the RL^Ga stoichiometry in 
their workii2iii*i£ Elemental analysis studies unavailable 
to these previous authors suggest that R2lr3Gag is the 
correct stoichiometry of these materials instead. The 
R2lr3Gag compounds with Ce and Yb show reduced mag- 
netic moments and the absence of magnetic order above 
0.04 Kil 

The thermodynamic properties of the single-impurity 
model have been calculated exactly using the Bethe- 
ansatz technique ^l^l^lZjl^ and also approximately 



within the Non-Crossing Approximation (NCA), which 
shows good agreement with the former. 19 However, to the 
best of our knowledge, it has always been assumed that 
the hybridization V m between any state of the magnetic 
configuration \m) and the conduction electrons is inde- 
pendent of m even when the CEF is included. 14 i 15 i 16 This 
is a requirement for the integrability of the problem^ al- 
though, there is no symmetry requirement to have the 
same V m for each orbital in the presence of a CEF. 
The different orientations of the orbitals relative to the 
nearest-neighbor atoms clearly indicate that the hy- 
bridizations must be a function of m. It is essential to 
consider this effect in order to explain the magnetic sus- 
ceptibility x(T) measurements of Yb2M 3 Gag (M= Rh, 
Ir) shown here as well as various other experimental ob- 
servations in related materials* 2 -!' 

Here, we report two examples of mixed-valence sys- 
tems, Yb2RIi3Gag and Yb2lr3Gag, in which the CEF and 
Kondo energy scales are of the same order of magnitude. 
We discuss a method based on the simple approximation 
scheme of Zwicknagl et al. for calculating dynamical and 
static properties for these types of systems^ The nov- 
elty of the method is that, in addition to the usual CEF 
effects, it also incorporates the important consequences 
of having different hybridization amplitudes V(m) be- 
tween the CEF orbitals and the conduction band. The 
calculated x{T) provides a quantitative description of the 
measured susceptibility in Yb2MsGag (M= Rh, Ir). We 
show that for Yb2h~3Gag a single hybridization V(m) is 
sufficient to describe the data, whereas for Yb2Rh 3 Gag, 
two hybridizations are necessary to adequately describe 
our results. 



I. EXPERIMENT 

Single crystalline rods of Yb2MaGag (T=Rh, Ir) with a 
tapered hexagonal morphology were grown by a Ga self- 
flux technique. Elemental analysis confirmed the correct 
2-3-9 stoichiometry. X-ray powder diffraction measure- 
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ments on crushed single crystals produce a spectrum that 
can be indexed in cither a hexagonal or C-centered or- 
thorhombic structure. Examination of the resulting lat- 
tice parameters suggests that the hexagonal unit cell is 
apparently a subcell of a larger orthorhombic cell (with 
&orthorhambic VS&hexagonai) ■ The room-temperature val- 
ues of the hexagonal lattice parameters are a — 7.471(6) 
A, c = 9.440(3) A for Yb 2 Rh 3 Ga 9 and a = 7.483(4) 
A, c=9. 441(2) A for Yb2lr3Gag, in agreement with 
Ref. |^ITlT|. The larger orthorhombic lattice constants 
are close to those reported for polycrystalline samples 
of the same compounds obtained by means of arc- or 
induction-meltingi 4 * 2 ^ 2 ^ X-ray single crystal diffraction 
studies (discussed in detail elsewherei£)that include the 
possible role of stacking fault disorder, as indicated in 
isostructural aluminum-rich compounds^, suggest that 
the samples crystallize in the hexagonal structure. Spe- 
cific heat and magnetization measurements were per- 
formed in a commercial (Quantum Design) PPMS and 
MPMS, respectively. The resistivity was measured using 
a standard four-probe technique, with the current paral- 
lel to the c axis in the temperature range of 0.5 — 300 K 
under zero applied field. 

Figure ^ shows the magnetic susceptibility, x(T), of 
Yb2Rli3Ga9 measured in H = 1 kOe. x(T) displays a 
broad maximum at 90 K (H || c) and 165 K (H _L c), typ- 
ical of mixed-valence compounds. The high-temperature 
magnetic susceptibility of Yb2Rh3Gag follows a Curie- 
Weiss law above 200 K yielding values of fi e ff = 4.3 
/is/Yb ; 9 = -50 K and n eff = 4.4 p B /Yb ; 6 = -293 
K for held parallel and perpendicular to c axis, respec- 
tively. The average effective magnetic moment is nearly 
the value of free Yb 3+ (fi e ff = 4.5/x.b). The marked dif- 
ference in the respective values for the Weiss temperature 
is due to the strong magnetocrystalline anisotropy. As 
will be discussed later, this anisotropy is also the origin 
of the orientation-dependent maxima in x(T). 

Figure [21 displays magnetic susceptibility data for 
Yb2lr3Gag. These data are consistent with stronger hy- 
bridization in Yb2li'3Gag as compared with Yb2Rh3Gag. 
The maxima in x(T) are shifted to higher temperature, 
and Curie- Weiss behavior is not fully recovered by our 
highest measurement temperature. As pointed out by 
Trovarelli et al.^ the upturns in susceptibility at lowest 
temperature may not be extrinsic. These authors reach 
this conclusion based on field-dependent magnetization 
measurements. That we see field-orientation dependent 
upturns at low temperature, especially in the case of 
Yb2Rli3Gag [x(T) for H || c is essentially temperature 
independent whereas x(T) for H _L c increases rapidly 
with decreasing temperature] adds further credibility to 
this assertion** 

The electrical resistivity p{T), measured with current 
applied along the hexagonal c-axis (the long axis of our 
rod- like crystals), of Yb2Rh-3Gag and Yb2lr3Gag is dis- 
played in the left insets of Fig. ^ and Fig. |3 respec- 
tively. The observed temperature dependences are char- 
acteristic of intermediate- valence metals and again reflect 



the higher characteristic temperature (i.e., the inflection 
point) in Yb2lr3Gag. Similarly, specific heat divided by 
temperature C/T versus T 2 is shown in the right in- 
sets of Figs, n and [21 f° r T=Rh, Ir, respectively. A 
low-temperature fit to C/T = 7 + /3T 2 gives values of 
7 = 45(25) mJ/mol-Yb K 2 and f3 = 0.75(0.91) mJ/mol- 
Yb K 4 , corresponding to a Debye temperature 9n — 250 
K (234 K) for M=Rh (Ir). The values of 7 reflect mod- 
erate mass enhancement in Yb2RIi3Gag and Yb2lr3Gag. 
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FIG. 1: Magnetic susceptibility x( T ) of Yb 2 Rh 3 Ga 9 . The 
lines are fits to the data using the model described in the text. 
Left inset: Electrical resistivity p(T) of Yb2Rli3Gag. Right 
inset: Specific heat divided by temperature C/T versus T 2 of 
Yb 2 Rh 3 Ga 9 . 



II. THEORY 

To account for the combined effects of Kondo hy- 
bridization and crystal electric field splitting, we use an 
approximation scheme developed in Refs. 



22 25] that we 

generalize to allow for different conduction electron hy- 
bridizations of each ground state multiplet of the mag- 
netic configuration. The approach is based on the non- 
crossing approximation (NCA) with the additional as- 
sumption that the density of states of the non-magnetic 
configuration by its pole at temperature T = 0: 



PoM = (1 - n f )S(uj - w ), 



(1) 



where 2 + n./ is the valence of Yb at T = 0, and uiq 
is the ground state energy, obtained by a variational 
ansatz that is exact for large degeneracy of the magnetic 
configuration^ The approximation Eq. 2] ceases to be 
valid at temperatures of the order of the charge transfer 
energy Eq — min(_E m ) + ep (see below). However, it has 
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FIG. 2: Magnetic susceptibility x(T) of Yb 2 Ir 3 Ga 9 . The lines 
are fits to the data using the model described in the text. 
Left inset: Electrical resistivity p(T) of Yb2rr3Gag. Right 
inset: Specific heat divided by temperature C/T versus T 2 of 
Yb2lr3Gag. 



led to very good agreement with results for the magnetic 
susceptibility of the full NCA and Bethe-ansatz in the 
isotropic casei^S Furthermore, it has the advantage that 
the results converge smoothly to those of the variational 
approximation for T = and is free from the usual arti- 
facts of the NCA at low 

The Hamiltonian is that of the impurity Anderson 
model including crystal and magnetic field: 



H = H Q + H 1 +H B + H h 



and 



H„ 



(2) 



with 



H 



£o|0)<0|, Hi =Y,Em\m)(m\, 

m 

km 

5^V T „(|m)<0|cfe TO + H.c.). 



(3) 



km 



Ho describes the ground state of the non-magnetic 4f 14 
configuration of Yb +2 ; Hi corresponds to the ground 
state multiplet of the 4f 13 configuration, distributed in 
four Kramers degenerate doublets belonging to three ir- 
reducible representations of the point group; c|, m creates 
a hole in an extended state with the same symmetry as 
the localized state \m). The density of these band states 
p is assumed constant and independent of to as usual. 
The term i? m i x allows for different hybridizations V m for 
each doublet. Hg describes the coupling of the magnetic 
configuration with an applied magnetic field B in the di- 
rection a. The value of g for Yb is 8/7. 



To be able to use the NCA at finite B when V m de- 
pends on to, we diagonalize first Hi + Hg and perform a 
canonical transformation on the ct„ in such a way that 
Hi + Hb takes the form of H\ with field dependent E m , 
and H m [ K retains the same form in the new basis. Calling 
\m(B)) the eigenstates of Hi + Hb, the transformation 
is: 



C ^ {B) = VjB) £^'K(0)||m(B))4 m ,, (4) 



where 



V>(B)=J2v*,\{m'(0)\\m(B))f 



(5) 



Working up to second order in B, diagonalizing first J a 
in each subspace of identical E m , denoting T m = irpV^, 
we obtain for E m (B) and T m (B) in the new basis: 



E m =E° m - gp B B(m\J a \m) + (gp B Bf £ 

v -r° +(ou m^ gV^aMf (6] 

lm-l m + (9^bB) (F o _ F0 \2 ' W 

where VL/ runs over all to' with E®, ^ E® and E®, = 

£ 4 Til til 1 '<v 

E m (0), = r m (0). In this way the Hamiltonian is 
mapped into H — Hb with field dependent parameters. 

Proceeding as in RefiSi, the ground state energy ojq(B) 
is obtained from the equation: 



t-^ \ ^ T m E^ 



ujo 



7T W + E n 



UJq 



(7) 



where we set the Fermi energy tp to zero and W is the 
width of the part of the conduction band occupied with 
holes. Using Xq(0) = —d 2 iv /dB 2 and neglecting 1/W in 
comparison with 1/Tq, where Tq = mhi(E m ) — luq is the 
stabilization energy of the correlated singlet, we obtain a 
closed expression for the susceptibility at T = 0: 



C«(0) = (5Ms) 2 £™/m[ 



x ( 



[m\J a \m) 
E? n ~ wo 



2^|(to|J q |to') 



El-Uo , ln fm-^0 }]) (g) 



El,-E« m (E° m -E°)* E° m ,-wo' 



where n./ TO , the occupation number of the state |to) at 
T = 0, is given by: 



Hfrn 



r„ 



n(E m - lu ) 



(9) 
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The same expression Eq. [5] is obtained finding first 
the variational wave function for B = and then using 
second order perturbation theory in B. While this pro- 
cedure is actually easier, it cannot be extended to T ^ 0. 

At finite T, the susceptibility is obtained from \a (T) = 
—d 2 F/dB 2 , where the free energy is given by: 



F = wo-Tln[(l-n/)(l+5^r T „7 m )], (10) 



with: 



I m = 



ir (u> + ljq - E m ) 2 + [r m (l - n f )f(-L 



(11) 
r(12) 



where f(uj) is the Fermi function. The second derivative 
of F is calculated using Eqs. l6l l$l and ITTI 

Up to now the theory corresponds to a single magnetic 
impurity coupled to band states. We find that in order 
to explain the observed magnetic susceptibility, in par- 
ticular the ratio of x(T = 0) and the maximum value of 
x{T), a small antiferromagnetic interaction between Yb 
ions should be included. In mean field, the susceptibility 
of the compound \% is given by: 



Xo 



Xo 



1 + IX* 



I = 



(.9MB ) 2 ' 



(13) 



where the sum runs over the exchange interactions of all 
sites that interact with the given one. 



III. DISCUSSION 

The best fits to the anisotropic magnetic susceptibility 
data for Yb2Rh3Gag and Yb2lraGag using our theoreti- 
cal framework are shown in Figs. ^and|21 respectively. 
The parameters associated with these fits are given in 
Table HJ To obtain these fits in practice, we use nr at 
T = B = rather than Eq and W as fitting parame- 
ters and determine ojq from Eqs. 10 and l|llfl . Because 
the upturns in susceptibility at low temperature are ei- 
ther extrinsic or periodic effects beyond the scope of our 
model, only data above 50 K were used in our fits. Sev- 
eral considerations influenced our fitting of the data in 
order to minimize the number of allowed free parameters. 
In order to obtain a larger susceptibility for B parallel to 
the c axis (x||)j states with larger angular momentum 
projection along that axis m c should lie lower in energy. 
However, the well-defined maximum in X|| a t intermedi- 
ate temperatures (T max ~100 K) points to an effective 
large degeneracy (this is clear from Bethe-ansatz results 
in the isotropic case^i) As a compromise, we took a four- 
fold degenerate ground state of H\ , containing the states 



with m c — ±7/2 and m c = ±5/2 along c (all belong- 
ing to the same irreducible representation). The remain- 
ing four states were also assumed degenerate. Allowing 
these quartets to split into closely spaced doublets does 
not qualitatively change the results. A larger admixture 
with the 4f 14 configuration (lower n/) decreases both sus- 
ceptibilities but increases relatively the upturn in \\\ as 
pointed out earlier^ 

The structure of the compounds suggest a larger hy- 
bridization for orbitals lying in the plane perpendicular to 
the c axis, which should correspond to higher m c . Thus, 
we begin our fitting procedure with two different values 
of r m , with the larger value corresponding to the ground 
state quartet. Increasing T m for both quartets leads to 
flatter curves and more similar values of x\\ an d X-L- As 
indicated in Table [fl a single hybridization parameter 

(i.e., r = r° ±7/2 = r° ±5 / 2 = r° 3/2 = r° 1/2 ) is sufficient 

to produce a high quality fit to the data for Yb2lr3Gag 
(Fig. [21). It is worth noting, however, that the included 
crystal field splitting is essential for describing the data. 

As the overall hybridization decreases from Yb2lr3Gag 
to Yb2Rli3Gag, the need for level-specific hybridization 
becomes apparent. In order to adequately fit the sus- 
ceptibility of Yb2Rli3Gag, especially the in-plane suscep- 
tibility (x±)j a quartet-specific hybridization had to be 
included. This is shown most clearly in Fig. [3] Here, we 
show best fits to the measured data using three calcu- 
lation approaches, that of the traditional Bethe-ansatz, 
the NCA with a single hybridization that was used for 
Yb2lr3Gag, and the NCA allowing for a different hy- 
bridization of the upper CEF quartet than that of the 
lower CEF quartet (T° 7 , 2 ±5 ^ 2 ^ r^_ 3 ^ 2 ±1 / 2 ). Our the- 
oretical approach also enables the calculation of 7, the 
low-temperature electronic contribution to specific heat. 
Using the best model parameters of Table we find 68 
mJ/mol Yb K 2 and 28 mJ/mol Yb K 2 for Yb 2 Rh 3 Ga9 
and Yb2lr3Gag, respectively, in reasonable agreement 
with measured values [45(25) mJ/mol-Yb K 2 for M=Rh 
(Ir)]. Inelastic neutron scattering measurements of the 
quasielastic linewidth and x-ray absorption spectroscopy 
measurements of n/ are presently in progress to further 
validate our models 



In summary, measurements of magnetic susceptibility, 
specific heat, and electrical resistivity have been per- 
formed on the mixed-valence compounds Yb2Rli3Gag 
and Yb2lr3Gag. Anderson impurity model calculations 
within the NCA approach describe the anisotropic mag- 
netic susceptibility indicating that it is essential to in- 
clude crystalline electric field effects. In Yb2lr3Gag, 
a single hybridization of the two split CEF quartets 
with the conduction electrons is needed to model the 
anisotropic x(T), while two different hybridizations of 
the two quartets are needed to fit the x(T) data of 
Yb 2 Rh 3 Ga9. 
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TABLE I: Fit parameters for the three calculations of the magnetic susceptibility of Yb2M3Gag (M=Rh, Ir) described in the 
text. Tk'- Kondo temperature [The values of the energy scale To for the NCA calculations have been scaled to compare with 



the Bethe-ansatz Kondo temperature Tk- The scale factor is given by: Tk = 



(2J+1) 
2irn f 



To, yielding the values of Tk listed in the 



table.]; I: molecular field constant; F: (single) conduction electron hybridization with CEF states; T^: conduction electron 
hybridization with CEF state |m); A: CEF splitting between upper and lower quartets; n/: /-occupation number at T = K; 

7: electronic specific heat coefficient; GOF for x II / -L c YliLi [ ^ X ' CT ^'*' ^ ] 2 /d, where N is the number of data points, 

a^ F is the error in x^ v , and d is the number of degrees of freedom!— The units of Tk, A, and T are [K]; units of / are [ m °^y b ]; 
units of 7 are f — , ™t „ ■ ! 1 ■ 

' L mol Yb Y- A 1 
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FIG. 3: (color online) Magnetic susceptibility x(T) of 
Yb2Rli3Ga9. The lines are fits to the data of the three calcu- 
lation approaches discussed in the text. 
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